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We have studied the pitch angles of spiral arms for 31 distant galaxies at z 0.7 from three Hubble 
Deep Fields (HDF-N, HDF-S, HUDF). Using the pitch angle - rotation velocity relation calibrated from 
nearby galaxies, we have estimated the rotation velocities of galaxies from the deep fields. These estimates 
have a low accuracy (~50 km/s), but they allow low-mass and giant galaxies to be distinguished. The 
TuUy-Fisher relation constructed using our velocity estimates shows satisfactory agreement with the 
actually observed relations for distant galaxies and provides evidence for the luminosity evolution of 
spiral galaxies. 



Q 1. Introduction 

^"-^ The pitch angle of spiral arms is a major parameter 
of the classical morphological classification of galaxies 
^^(Hubble 1936). The pitch angle is the angle between 
O the tangents to the spiral arm and to the circumfer- 
ence centered at the galaxy nucleus drawn through a 
given point. Galaxies with tightly wound spirals and 
open arms have small and large pitch angles, respec- 
1-H tively. Typically, the pitch angle lies within the range 
> r^Q° - 30° (Kennicutt 1981; Ma 2001). 
I , In 1981, Kennicutt published a paper in which the 
relation between the pitch angle of a spiral galaxy 
. and its maximum rotation velocity was studied. It 
followed from this paper that galaxies with tightly 
^—1 ■(vound arms (smaller pitch angles) rotated, on aver- 
^""^ age, faster than those with more open arms (larger 
J> pitch angles), with this dependence being linear. In 
^> recent years, investigating the pitch angle of spiral 
^ arms has become increasingly popular, because sev- 
. 5^ eral unexplained empirical relationships of this quan- 
tity to the parameters of the galaxy rotation curve 
and to the mass of the central black hole have been 
found (see, e.g., Seigar et al. 2008; Shields et al. 2010). 

The goal of this paper is to develop the results from 
Kennicutt (1981) and to try to apply the pitch angle 
- rotation velocity relation to determine the rotation 
velocities of distant spiral galaxies from the Hubble 
Deep Fields HUDF, HDF-N, and HDF-S. Since many 
of these galaxies are very faint and have small angular 
sizes, measuring their rotation velocities by spectro- 



scopic methods requires instrumentation that will not 
be available in the near future. However, knowledge 
of these velocities is needed to study the evolution 
of galaxies at high z. In addition, there are virtually 
no data on the pitch angles of spiral arms for distant 
spiral galaxies at present. 

All numerical values in the paper are given for the 
cosmological model with the Hubble constant 70 km 
s-i Mpc-\ = 0.27, and JIa = 0.73. 

2. Measurements 

2.1. Determining the orientation of galaxy disks 

An important preliminary stage of our work to study 
the shape of spiral arms is to estimate the orientation 
of galaxy disks in space. This orientation is speci- 
fied by two parameters: the inclination i of the disk 
plane to the plane of the sky and the position angle 
PA of the major axis. Knowing these parameters is 
needed for two reasons. First, because of the galaxy 
inclination to the plane of the sky, the apparent spi- 
ral structure is distorted and, therefore, the galaxy 
image should be deprojected to the "face-on" orien- 
tation. Second, the galaxy rotation velocity derived 
from spectroscopic measurements should also be cor- 
rected for the inclination. 

In this paper, we decided to use the relatively new 
method of spiral-arm monotony (SAM; Poltorak and 
Fridman 2007; Fridman and Poltorak 2010). This 
method is based on the assumption that every spiral 
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arm is a monotonic function, i.e., following along the 
spiral from the galaxy center to the periphery, the 
radius must increase monotonically > 0, where 
r is the distance from a point in the spiral arm to 
the galaxy center and (j) is the azimuthal angle). In 
this case, the spiral projected onto the plane of the 
sky can be represented by a nonmonotonic function. 
Thus, the domain of i and PA for which the depro- 
jected spiral is monotonic will be the domain of pos- 
sible inclinations and position angles. 

An example of using the SAM method to estimate 
the orientation of four nearby galaxies is shown in 
Fig. 1. We see from the figure that the domains of 
possible values are relatively small, which allows both 
inclination and position angle of the galaxy major 
axis to be estimated with a good accuracy. 

2.2. The pitch angle 

Once the galaxy images have been corrected for the 
inclination, we can turn to the determination of 
spiral-arm pitch angles. As the published data show, 
the pitch angles for the same objects often differ 
markedly. Therefore, for reliability, we decided to im- 
plement two different, completely independent meth- 
ods and to compare their results. 

The first method (below called the interactive one) 
is based directly on the search for the angle between 
the tangents to the spiral arm and to the circumfer- 
ence centered at the galaxy nucleus. If the spiral arm 
is represented in polar coordinates: r = r{(j)) (the co- 
ordinate origin at the galaxy center), then the pitch 
angle for a logarithmic spiral can be determined from 
the formula (see, e.g., Binney and Tremaine 1987) 









dr 



(1) 



The main steps in determining the pitch angle are: 
(1) determining the coordinates of the galaxy cen- 
ter; (2) finding the polar coordinates of several (~10) 
points on the arm; (3) applying pairwise Eq. (1) to the 
derived coordinates to obtain the set of pitch angles 
corresponding to different arm segments; (4) averag- 
ing the pitch angles to obtain the final result. 

To reduce the random errors, this procedure was 
repeated several times for each arm. If the galaxy 
had several (as a rule, two) large-scale arms suitable 
for measurements, then their pitch angles were also 
estimated in the same way. As the final pitch angle 
for a given galaxy, we took the value averaged over 
the measured arms. 



The second method is based on a Fourier analysis 
of the distribution of points in the spiral arms of a 
galaxy (Considere and Athanassoula 1982). If the dis- 
tribution of points in the galaxy arms is represented 
as the sum of delta functions of their polar coordi- 
nates: 



1 ^ 

— ^5(it - Ui)5{(p - (pi), 
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where Ui = In(rj), then its Fourier transform 
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gives the coefficients in the expansion of this distribu- 
tion in terms of logarithmic spirals. The pitch angle 
can be found from the formula 



^ = arctan 



m 



(4) 



where m is the number of spiral arms in the galaxy 
and Pmax is the value of the parameter p at which the 
function |A(p, m)| has a maximum. 

3. Measuring the pitch angles for galaxies 
of a local sample 

3.1. The galaxy sample 

To study the spiral pattern of nearby galaxies, 
we used the galaxy sample described by Kennicutt 
(1981). The sample includes 113 spiral galaxies of var- 
ious types with a clearly distinguishable spiral struc- 
ture. Since the arms of barred galaxies are described 
more poorly by a logarithmic spiral, most of the sam- 
ple galaxies are without bars or have small bars. 

The galaxy images in the fits format were retrieved 
from the NEE0, which contains references to observa- 
tions with different instruments and in different spec- 
tral ranges. The pitch angle depends on the filter; for 
our purposes, we used only the i?-band images, be- 
cause the spiral structure is seen best in blue bands. 
No photometric calibration of the galaxy images was 
required and, therefore, we used the data obtained 
with different instruments; from several possible op- 
tions, we chose the best one from the viewpoint of 
visual image reduction (where possible, a higher res- 
olution and a better view of the spiral structure). 
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Fig. 1. Results of the operation of the spiral-arm monotony method for the galaxies NGC 2336, NGC 3294, 
NGC 4939, and NGC 6118 as an example. The black color indicates the domain of possible inclinations and 
position angles. 



The maximum rotation velocities of the sample 
galaxies Vmax found from the width of the HI (A=21 
cm) line profile were retrieved from the HyperLEDAl. 
The HyperLEDA rotation velocity have already been 
corrected for the inclination found from the appar- 
ent flattening of galaxies and, therefore, we initially 
eliminated this correction and subsequently applied it 
again using the galaxy inclination estimated by the 
SAM method (see below). 

3.2. Measuring the pitch angle 

For all galaxies from the sample by Kennicutt (1981), 
we estimated their orientation parameters, the in- 
clination and position angle of the major axis, by 
the SAM method. The agreement of our inclina- 
tions with those in HyperLEDA is, on average, good: 
{iMCP — ilEDA) = —2° lb 13°. However, for sev- 
eral galaxies, the apparent isophote flattening method 
gives, obviously, incorrect values due to various kinds 
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of peculiarities in the shapes of the outer isophotes of 
galaxies. 

The SAM method has a limitation at low inclina- 
tions (i.e., when the galaxies are seen almost face-on): 
the error in the i and PA estimates increases greatly 
with decreasing inclination. For this reason, we failed 
to estimate the orientation parameters for 16 galaxies 
and they were excluded from our sample, in which, 
thus, 97 objects remained. 

Once the images have been deprojected, we deter- 
mined the spiral-arm pitch angle for each galaxy by 
the two methods. At this step, we excluded several 
more galaxies with an excessively irregular and asym- 
metric structure from our sample. We failed to find 
the maximum rotation velocities for four galaxies. In 
addition, we decided to restrict ourselves to the galax- 
ies with inclinations in the range 30° < i < 60°, be- 
cause both pitch angles and maximum rotation ve- 
locities can be reliably estimated for such galaxies 
and their pitch angle - rotation velocity relation (see 
below) is seen much better. As a result, 46 galaxies 
remained in our sample of nearby galaxies with pitch 
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angles measured by the two methods and with known 
rotation velocities, 43 and 3 of which have two-armed 
and three-armed spiral patterns, respectively. 

Figure 2 compares our pitch-angle measurements 
with the results from Kennicutt (1981) for both meth- 
ods. The straight lines indicate linear fits to the data. 
They suggest the existence of certain systematics in 
the pitch-angle measurements probably attributable 
to the difference between the methods used. The 
mean differences between the measurements for the 
same galaxies are small: A// = 0.2° zt 3.7° (the dif- 
ference between the pitch angle measured by the 
first (interactive) method and the angle measured by 
Kennicutt) and A/x = 0.0° it 3.5° (the same for the 
second method). 

Comparison of the two pitch-angle determination 
methods described above shows good mutual agree- 
ment between the results: the difference between the 
angles estimated by the Fourier and interactive meth- 
ods is Afi = —0.2° lb 3.0°. Both methods yield sim- 
ilar results, although, formally, the accuracy of the 
method based on a Fourier analysis is higher. 

Figure 3 shows the logarithmic spirals constructed 
from the spiral-pattern parameters we determined 
and, given with the galaxy inclinations and position 
angle, they were superimposed on the observed im- 
ages of NGC 2997 and NGC 4254. 

3.3. The pitch angle - rotation velocity relation 

for nearby galaxies 

Figure 4 shows the pitch angle - rotation velocity 
relations constructed by the two methods. As we 
see from the figure, both approaches give significant 
correlations. However, the relation constructed from 
the measurements using the Fourier transform has a 
higher statistical significance (the linear correlation 
coefficient for it is 0.725 versus 0.600 for the relation 
constructed by the interactive method). In our sub- 
sequent discussion, we will use the pitch angles found 
by a Fourier analysis. Note also that the galaxies with 
a three-armed pattern in Fig. 4 are near the average 
relation for galaxies with a two-armed pattern. 

The final empirical relation found from 46 nearby 
galaxies (the solid straight line in Fig. 4b) is 

/x(°) = (-0.049 ± 0.008) • y^„^(km/s) (5) 

+ (22.85 ±1.76). 

Curiously, our Galaxy may satisfy this relation. It 
follows from Eq. (5) that /x = 12° for the Milky Way 



with Vmax = 220 km/s. This pitch angle is close to 
the estimates of various authors (see, e.g. Tables 1 
and 2 in Vallee 2005). However, greatly differing pitch 
angles of the Milky Way spiral pattern are also en- 
countered in the literature (see, e.g., Levine et al. 
2006). 

4. Measuring the pitch angles for galaxies 
from the deep fields 

4.1- The galaxy sample 

To find distant galaxies with a spiral structure, we 
used the original frames of the Hubble Deep Fields 
HDF-N and HDF-S (Ferguson et al. 2000) as well as 
the Ultra Deep Field HUDF (Beckwith et al. 2006). 
The main problem in compiling the sample was that 
the angular sizes of many distant galaxies are too 
small and, therefore, several objects with a distinct 
spiral structure were not included in the sample to 
avoid the errors due to image discreteness. In addi- 
tion, a considerable number of deep-field galaxies ap- 
pear asymmetric and peculiar and no regular spiral 
arms can be drawn for them. 

Our final sample includes 31 galaxies (20 in HUDF, 
six in HDF-N, and five in HDF-S) with a clearly dis- 
tinguishable spiral pattern, rule, a two-armed 
ones as in the case of nearby galaxies. The mean red- 
shift of these galaxies is (z) = 0.69 ± 0.30. Our list of 
distant galaxies is presented in the table. 

4-2. Measurements 

For all sample galaxies, we downloaded their images 
in the F606W filter, because this is the only com- 
mon filter for all three fields. To estimate the orien- 
tation of distant galaxies, we used the SAM method. 
In the cases where this method jdclded no definite 
results, we estimated the inclination and orientation 
from the outer isophotes of the galaxies. (Such ob- 
jects were discarded among the nearby galaxies but 
were retained among the distant galaxies, because of 
the small size of their sample. Their inclinations are 
marked by the asterisk in the table.) 

The results of our measurements of the inclination 
(i) and spiral-pattern pitch angle (/i) are summarized 
in the table. If a galaxy has only one arm suitable for 
measurements, the results are presented only for it. 
The galaxy redshifts (the fourth column of the table) 
were taken from Wolf et al. (2004) (the COMBO-17 
project) for HUDF, from Wirth et al. (2004) for HDF- 
N, and from Sawicki and Mallen-Ornelas (2003) for 
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Fig. 2. Comparison of our pitch angles with those from Kennicutt (1981): (a) for the interactive method; 
(b) for the method based on a Fourier analysis of galaxy images. 




Fig. 3. Images of the galaxies NGC 2997 (left) and NGC 4254 (right) with the logarithmic spirals corre- 
sponding to the pitch angles measured by the interactive method superimposed on them. 



HDF-S. (We used the redshift for the HUDF galaxy 
No. 16 from Coe et al. (2006), because z from Wolf 
et al. (2004) leads to an unrealistically low luminos- 
ity of the galaxy.) The eighth column of the table 
presents the maximum rotation velocity found from 
the empirical relation (5) and its error. We estimated 
the rotation-velocity error using the formula for the 
propagation of the mean error (the formula for the er- 
ror of a function of several variables) from the pitch 



angle measurement error and the errors of the numer- 
ical coefficients in Eq. (5). The fifth column gives the 
absolute magnitudes of the galaxies in the Hubble 
Space Telescope F606W filter. To find the absolute 
magnitudes, we used the ^-corrections for galaxies of 
the corresponding types from Bicker et al. (2003). 
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Fig. 4. Pitch angle - maximum rotation velocity relations derived by the two methods: the interactive method 
(a) and Fourier analysis (b). The solid lines are linear regressions; the dashed lines are iter deviations. 



5. Comparison of the local and distant 
galaxies 

Figure 5 presents the histograms showing the pitch- 
angle distributions for the local-sample galaxies and 
for the galaxies from the deep fields. We see good 
agreement in the distribution of angles: a maximum 
near 14° and a gradual decline on both sides of the 
maximum. Thus, one might expect the distributions 
of rotation velocities for galaxies at different z to be 
similar. 

Using the maximum rotation velocities of distant 
galaxies estimated from the empirical relation (5) 
(the table), we constructed their Tully-Fisher (TF) 
relation (Fig. 6) . The solid line in the figure indicates 
the relation for nearby galaxies in the same color band 
(Sakai et al. 2000). As wc see from the figure, the dis- 
tant galaxies are located on this plane with approxi- 
mately the same slope as that for the nearby galaxies, 
but their distribution appears shifted towards higher 
luminosities. The value of this shift depends on the 
galaxy redshift: the galaxies at z < 0.66, on average, 
follow the local relation, while those at z > 0.66 are 
shifted upward (Fig. 6). 

The most probable cause of the systematic shift in 
the positions of distant galaxies is the evolution of 
their luminosity. As was shown in numerous works 
devoted to studying the evolution of the TF relation 



Local galaxies 



Distant galaxies 



Pitch angle (degrees) 



Fig. 5. Histograms representing the pitch-angle dis- 
tributions of galaxies: for nearby galaxies (top) and 
for distant galaxies from the deep fields (bottom). 
The fractions of galaxies in the corresponding bins 
are shown along the vertical axes. 



and the galaxy luminosity function, spiral galaxies at 
z ^ 1 are brighter than nearby galaxies with the same 
maximum rotation velocity by ~ 1"* (see, e.g.. Table 1 
in Portinari and Sommcr-Larscn 2007). The right 
panel of Fig. 6 shows the TF relation corrected for 
the luminosity evolution of distant galaxies (the evo- 
lution correction was applied as prescribed by Bicker 
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Table 1. Galaxies from the deep fields 
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et ah 2003). As we see from the figure, allowance for 
the luminosity evolution slightly improved the agree- 
ment between the relations for distant and nearby 
galaxies. 

Consider two subsamples approximately equal in 
size - galaxies with z < 0.66 (15 objects) and z > 0.66 
(16 galaxies). The mean characteristics of the first 
subsample are {z) = 0.44±0.04 (the standard error of 
the mean is given), {MpeoGw) = —20.91 =b 0.36 (only 
the fc-correction was applied), and {Vmax) = 187 ±20 
km/s; the characteristics of the second subsample are 
{z) = 0.93 ± 0.05, {MFdoew) = -21.87 ± 0.22, and 
{Vmax) = 167 ± 16 km/s. 



Using the local TF relation from Sakai et al. (2000), 
we can estimate the expected absolute magnitude for 
a galaxy with Vmax = 187ib20 km/s to be Mpgoew = 
—20.7 lb 1.0 and for a galaxy with Vmax = 167 it 16 
km/s to be Mpqqqw = —20.3 it 0.9. The mean ob- 
served magnitude for galaxies with (z) = 0.44 is 
brighter than the expected one by AM = 0."*2=b 1.^*1 
and for the more distant subsample by AM = l."^6ib 
0.™9. According to the model by Bicker et al. (2003), 
the corresponding values of the luminosity evolu- 
tion for Sb-Sc galaxies are 0."*3-0.'"5 and 0.'"6-0."*9. 
Given the approximate nature of our method for esti- 
mating the rotation velocities and the small size of the 
galaxy sample, it may be concluded that the shift of 
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Fig. 6. The Tully-Fisher relation for distant galaxies derived using the pitch angle - rotation velocity relation: 
the absolute magnitudes were calculated by applying only the ^-correction (left) and the correction for 
evolution (right). The straight lines correspond to the local relation (Sakai et al. 2000); the dashed straight 

lines indicate the ztSo" deviations. The circles, triangles, and squares represent the galaxies from HDF-N, 
HDF-S, and HUDF, respectively. The open and filled symbols indicate the galaxies with z < 0.66 and 
z > 0.66, respectively. 



distant galaxies in the TF relation observed in Fig. 6 
agrees satisfactorily with that expected from the lu- 
minosity evolution with z. 

6. Conclusions 

We determined the spiral-arm pitch angles for 46 
nearby galaxies with 30° < i < 60° by two differ- 
ent methods (the interactive one and using a Fourier 
analysis of images). The inclination of the galaxy 
plane to the line of sight was found by the relatively 
new method of spiral-arm monotony (Poltorak and 
Fridman 2007; Fridman and Poltorak 2010). 

We confirmed the conclusion by Kennicutt (1981) 
about the existence of a significant correlation be- 
tween the rotation velocity of a galaxy and the pitch 
angle of its spiral arms. This correlation is best traced 
for intermediate galaxy disk inclinations, when both 
pitch angles and rotation velocities can be found with 
a good accuracy. 

We measured the pitch angles for 31 spiral galax- 
ies from several Hubble Deep Fields at mean redshift 
(z) ~ 0.7. Using the local empirical pitch angle - rota- 
tion velocity relation, we estimated the maximum ro- 
tation velocities of distant galaxies (the table) . These 



estimates have a low accuracy (Ri 50 km/s, see the 
table) , but they allow low- mass and giant galaxies to 
be confidently distinguished. 

We constructed the TF relation between the ab- 
solute magnitudes and maximum rotation velocities 
of distant galaxies estimated from the pitch angles of 
the spiral pattern. Despite the large scatter of data, 
we can tentatively conclude that the distant galaxies 
follow the local TF relation with approximately the 
same slope and, in addition, show evidence of lumi- 
nosity evolution. These results agree with the direct 
measurements of the TF relation for distant galaxies 
by spectroscopic methods. 

The "morphological" estimates of the rotation ve- 
locities for distant galaxies can be useful in study- 
ing the evolution of spiral galaxies, especially in the 
cases where no data on the rotation of a galaxy can 
be obtained by spectroscopic methods because of its 
faintness or its visibility conditions (the galaxy is seen 
almost face-on). 
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